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ABSTRACT 

We present time-resolved and phase-resolved variability studies of an extensive X-ray high- 
resolution spectral dataset of the 6 Ori Aa binary system. The four observations, obtained with 
Chandra ACIS HETGS, have a total exposure time of f«479 ks and provide nearly complete 
binary phase coverage. Vaiiability of the total X-ray flux in the range 5-25 A is confirmed, with 
maximum amplitude of about ±15% within a single «125 ks observation. Periods of 4.76d 
and 2.04d are found in the total X-ray flux, as well as an apparent overall increase in flux level 
throughout the 9-day observational campaign. Using 40 ks contiguous spectra derived from the 
original observations, we investigate variability of emission line parameters and ratios. Several 
emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, 
variations of the X-ray emission line widths as a function of the binary phase are found in a 
binary system, with the smallest widths at (/)=0.0 when the secondary 5 Ori Aa2 is at inferior 
conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission 
line width variability to the presence of a wind cavity created by a wind-wind collision, which 
is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary 
wind, thus producing phase locked X-ray variability. 

Subject headings: start: individual (i5 Ori A)—binaries: close — binaries: eclipsing — stars: early-type 
— stars: X-ray spectroscopy—stars: X-ray variability 
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1. INTRODUCTION 


Stellar winds of hot massive stars, primar¬ 
ily those with M > 8 Mq, have important ef¬ 
fects on stellar and galactic evolution. These 
winds provide enrichment to the local interstel¬ 
lar medium via stellar mass loss. On a larger 
scale, the cumulative emichment and energy 
from the collective winds of massive stars in 


a galaxy are expecfed fo play a pivofal role 

in driving galactic winds ( 

^eifherer el al. 1992: 

Oooenheimer & Dave 

2006 


McKee & Osfriker 


20071). The number of massive stars in any star- 


forming galaxy, as well as their tendency to be 
found in clusters, are critical parameters for de- 
termining a gala xy’s energy budge t and e volution. 


Oskinoval (120051) and lAgeitz et al.l (l2013h showed 


that the energy from the winds of massive stars 
will dominate over the energy from supernovae in 
the early years of massive star cluster evolution. 
While substantial progress has been made over 
the last several decades in modeling massive star 
winds (IPuls et al.lll996h . many questions remain, 
such as the degree of clumping of the winds, the 
radial location of different ions and temperature 
regimes in the wind with respect to the stellar sur¬ 
face, and the origin of Corotating Interaction Re¬ 
gions (CIRs) representing large-scale wind pertur¬ 
bations. X-ray observations have provided power¬ 
ful diagnostic tools for testing models, but a fully 
consistent description of the detailed structure of 
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a stellar wind is still elusive. 

Variability in the winds of massive stars can 
be an important probe of the structure of the stel¬ 
lar winds. There can be multiple causes of X- 
ray variability in massive stars. Large-scale struc¬ 
tures in the winds, as trac ed by Discrete Ab sorp- 
tion Components (DACs, iKaper et all 119990 and 
possibly linked to CIRs, may be associated with 
shocks in the wind and thereby potentially af¬ 
fect the X-ray emission. X-ray variations of this 


Pup, and ^ Per (Oskinova el al. 

2001; 

Naze ef al. 

2013; 

Massa el al. 

2014 

). Also, X-ray variabilily 


with the same period as but larger amplitude than 
known pulsational activity in the visible domain 


was r ecently detected in CMa (lOskinova et al. 


20141) and possibl y in the hard band of /? Cru 


(ICohen et al.ll2008h . The exact mechanism giving 
rise to these changes remains unclear. Notably, 
other pulsating massive stars do not show such 


X-ray “pulsations,” such as 0 Cen (IRaassen et al. 


20051) and /? Cep (IFavata et al.l 120091) . Smailer- 


scale structures such as clumps can also produce 
X-rays, albeit at lower energies than large-scale 
structures. It is also possible that some X-ray 
variations in massive stars are stochastic in nature 
and are not correlated with any currently known 
timescale. 

Another cause for X-ray variability is possible 
in magnetic stars. When a strong global magnetic 
field exists, the stellar wind is forced to follow the 
field lines, and fhe wind flowing from fhe fwo stel¬ 
lar hemispheres ma y then collide at fhe e quator , 
generating X-rays (IBabel & Monfmerlel 119970 . 
Such recurrenf variations have b e en defecfed 
in 0^ Ori C ( Sfelzer ef al. 2005Gagne ef al. 


20051) . HD 191 612 (INaze ef al.ll2010h . and pos¬ 


sibly Trl6-22 (INaze ef all 120141) . alfhough the 
absence of large variations in the X-ray emis¬ 
sion of the magnetic star r Sc o is puzzling in 
this context dlgnace et al.l 120101) . Even stronger 
but very localized magnetic fields could also be 
presenf, e.g. associafed wifh brighf spofs on fhe 
sfellar surface fhaf are required fo create CIRs 
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dCranmer & Qwockill 19961 : ICantiello et alJl2009h . 

In multiple systems, the collision of the wind 
of one star with th e wind of ano t her ca n produce 
X-ray variations dStevens et alJ Il992h . Wind- 
wind collision emission may vary with binary 
phase, with inverse distance in eccentric sys¬ 
tems, or due to changes in line o f sight ab¬ 


sorpti on, as observ ed in HD 93403 dRauw et al. 


2002h. Cvg OB2 9 dNaze et alJl2012r) . V444 Cyg 


( Lomax et alJ l2015h. and possibly HD 93205 
( Antokhin et al. 2003 1. 

Delta Ori A (Mintaka, HD 36486, 34 Ori) is 
a nearby multiple system that include s a close 
eclips ing binary, S Ori Aal (09.5 II: IWalborn 


I972h and 6 Ori Aa2 ((BI V: IShenar et aLll20I5h. 


herein Paper IV), with pe riod «5.73d dHarvin et al. 


20021: iMaver et al.ll2010l) . This close binary is or¬ 


bited by a more distant com panion star, 6 Ori Ah 
(f«B0 IV: Pablo et al. ( 2015 ). herein Paper III; Pa¬ 
per IV), having a period of ^^346 yr ( Heintzlll987 ; 
Perryman & ESAI 1997 : Tokovinin 2014 ). The 


components Aal and Aa2 are separated by about 
43 R0(2.6 RAai, Paper III), and the inclination 
of Ri76°±4° (Paper III) ensures eclipses. We ac¬ 
quired ?«479 ks of high resolution X-ray grating 
spectra with a Chandra Large Pro gram to observe 


nearly a full period of 5 Ori Aa dCorcoran et al. 


(120151), herein Paper I). Simultaneous with the 
acquisition of the Chandra data. Microvariability 
and Oscillations of STars (MOST) space-based 
photometry and ground-based spectroscopy at nu¬ 
merous geographical locations were obtained and 
are reported in Paper III. 

Previous X-ray observations of the 6 Ori A sys¬ 
tem from Einstein show ed no significant variabil¬ 
ity ( Grady et al. 1984). ROSAT data for <5 Ori 
A were studied by iHaberl & Whitel d 19931) . who 
found modest 2 - 0 - variability but no obvious phase 
dependence; the ICorcoranI d 19961) reanalysis of 
the ROSAT data showed similar results. A single 
previous Chandra HETGS obseryation of 6 Ori 
Aa was analyzed by iMiller et all d2002l) . Fitting 
the emission lines using Gaussian profiles, they 
found the profiles to be symmetrical and of low 


FWHM, considering the estimated wind veloc¬ 
ity. The 60 ks expos ure time covers abo u t 12% 
of the orbital period. iRaassen & PollockI d2013h 
analyzed a Chandra FFTGS observation with ex¬ 
posure time of 96 ks, finding some variability in 
the zeroth order image; they were not able to de¬ 
tect any variability in the emission lines between 
two time splits of the observation. 


This paper is part of a series of papers. Other 
papers in this series address the parameters of the 
composite Chandra ^479 ks spectrum (Paper I), 
the simultaneous MOST and spectroscopic obser¬ 
vations (Paper III), and UV-optical-X-ray wind 
modeling (Paper IV). In this paper (Paper II), 
we investigate variability in the X-ray flux in the 
Chandra spectra. Sect. |2] describes the Chan¬ 
dra data and processing techniques. Sect. |3] dis¬ 
cusses the overall X-ray flux variability of the ob¬ 
servations and period search. Time-resolved and 
phase-resolved analyses of emission lines are pre¬ 
sented in Sect. |4l In Sect. [5] we relate our results 
of phase-based variable emission line widths to 
a colliding wind model developed for this binary 
system in Paper I, and discuss possible additional 
sources of variability in 6 Ori Aa. Sect, [^presents 
our conclusions. 


2. OBSERVATIONS AND DATA REDUC¬ 
TION 

Delta Ori Aa was observed with the Chandra 
Advanced CCD Imaging Spectrometer (ACIS) 
instrument using the High Ener gy Transmission 


Grati ng Spectrometer (HETGS) (ICanizares et al. 


2005h for a total exposure time of ?«479 ks, cov¬ 
ering parts of 3 binary periods (see Table [2] for a 
list of observations and binary phases). Four sep¬ 
arate observations were obtained within a 9-day 
interval. The HETGS consists of 2 sets of grat¬ 
ings: the Medium Energy Grating (MEG) with 
a range of 2.5-31 A (0.4-5.0 keV) and resolution 
of 0.023 A FWHM, and the High Energy Grat¬ 
ing (HEG) with a range of 1.2-15 A (0.8-10 keV) 
and resolution of 0.012 A FWHM. The resolu¬ 
tion is approximately independent of wavelength. 


3 





























































































Table 1 

System parameters for S Ori Aa1+Aa2. 


Parameter 

Value 


Sp. Type (Aal) 

09.511 

Sp. Type (Aa2) 


D [pc] 

380 (adopted)® 

RIRq] (Aal) 

16.5 ± 1 

R[Rq] (Aa2) 

6-5l?,5 

Binary Period “ 

P[d] 

5.732436'' 

Eo (primary min, HJD) 

2456277.790 ± 0.024 

Tq (periastron, HJD) 

2456295.674 ± 0.062 

“[-R©] 

43.1 ± 1.7 

i [deg.] 

76.5 ±0.2 

Ld [deg.] 

141.3 ±0.2 

id [deg. yr~^] 

1.45 ± 0.04 

e 

0.1133 ±0.0003 

7 [km s~^] 

15.5 ±0.7 

Periastron-based 4> 

0.116-i-Photometric-based 4> 

MOST optical secondary periods [d] 

MOSTpi 

2.49 ±0.332 

MOSTp2 

4.614 ± 1.284 

MOSTps 

1.085 ±0.059 

MOSTp4 

6.446 ±2.817 

MOSTps 

3.023 ± 0.503 

MOSTpe’^ 

29.221 ± 106.396 

MOSTpr 

3.535 ±0.707 

MOSTps 

1.01 ±0.051 

MOSTpg 

1.775 ±0.162 

MOSTpio 

2.138 ±0.24 

MOSTpii 

1.611 ±0.133 

MOSTp42 

0.809 ± 0.032 

MOSTpis 

0.748 ± 0.027 


^Shenar et al. (2015) 

^from the low-mass model solution of Pablo et al. (2015) 
iMalz ADellaniz~^ ail <201. 

"^Mayer et al. (2010) 

®This peak is likely an artifact due to a trend in the data. It is not consid¬ 
ered real, but it is formally significant and included in the fit. 
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The Chandra ACIS detectors record both MEG 
and HEG dispersed grating spectra as well as the 
zeroth order image. Due to spacecraft power 
considerations, it was necessary to use only 5 
ACIS CCD chips instead of the requested 6 for 
these observations. Chip S5 was not used, mean¬ 
ing wavelengths longer than about 19 A in the 
MEG dispersed spectra and about 9.5 A in the 
HEG dispersed spectra were only recorded for the 
“plus” side of the dispersed spectra, reducing the 
number of counts and effective exposure in these 
wavelength regions. The standard data products 
distributed by the Chandra X-ray Observator 
were further processed w ith TGCat softwan 
(IHuenemoerder et al.l l201lh . Specifically, each 
level 1 event file was processed info a new level 
2 evenf file us ing a package of sfan dard CIAO 
analysis fools ( Eruscione et al. 2006h . Addition¬ 
ally, appropriate redistribution matrix (RME) and 
area auxiliary response (ARE) files were calcu- 
lafed for each order of each specfrum. TGCaf 
processing produced analysis producfs wifh sup¬ 
plemental statistical information, such as broad- 
and narrow-band count rates. 


MOST photometry observations were obtained 
of 6 Ori Aa for approximately 3 weeks, including 
the 9 days of the Chandra observations. Eig. [J 
shows the simultaneous MOST lightcurve aligned 
in time to the Chandra lightcurve. The Chandra 
lightcurve in this figure is fhe ±1 orders of fhe 
HEG and MEG combined in fhe 1.7 < A < 25.0 
A range, binned af 4 ks, wifh Poisson errors. Eig. 
|2] shows fhe same dafa ploffed wifh binary phase 
rafher fhan time. The MOST lighfcurves from 
several orbifs are averaged and overploffed in fhe 
figure fo show fhe optical variabilify. 

Table [Ulisfs fhe specfral fypes and radii of fhe 
(5 Ori Aal and Aa2, as well as fhe orbifal parame- 
fers and fhe MOST secondary periods. In fhis pa¬ 
per 0=0.0 refers fo fhe binary orbifal period and 
denofes fhe time when fhe secondary is in fronf 
of fhe primary (deeper optical minimum) and 0 = 


^available for public download: http://tgcat.mit.edu 


0.5 denofes fhe time when fhe primary is approx- 
imafely in fronf of fhe secondary (shallower opti¬ 
cal minimum). While primary minimum is a def¬ 
inition, fhe secondary minimum af 0=0.5 is only 
approximafe, since the orbit is slightly elliptical 
and also varies slowly with apsidal motion. The 
actual secondary minimum is currently 0=0.45 
(Paper III). Actual current quadrature phases are 
0=0.23 and 0.78. To avoid confusion, we use the 
phases in this paper that would be assumed with 
a circular orbit (i.e., 0=0 for inferior conjunction, 
0=0.5 for superior conjunction, 0=0.25 and 0.75 
for quadrature). A lso, we use the ephemeris of 


Maver et all (120101) . No evidence of X-ray emis¬ 


sion from the tertiary star was seen in any of the 
Chandra observations (Paper I), so the spectra 
represent only 6 Ori Aal and Aa2. 


3. OVERALL VARIABILITY OF X-RAY 
FLUX 


The lightcurve of the dispersed Chandra spec¬ 
trum of 6 Ori Aa, shown in lower part of Fig. [J 
shows the spectrum was variable throughout in 
X-ray flux wifh a maximum amplifude of abouf 
±15% in a single observation. During fhe firsf and 
second observations fhe X-ray flux varied by Ri5- 
10%, followed by a ^15% decrease in fhe fhird 
and a R:^15% increase in fhe fourth. Nofe fhaf none 
of fhe X-ray minima in fhe residual lighfcurve 
aligns wifh an opfical eclipse of fhe 6 Ori Aa sys¬ 
tem. Fig. [J suggests an increase in overall X-ray 
flux wifh time. From fhe beginning fo end of fhe 
9-day campaign fhere is a R:^25% increase in fhe 
mean counf rate. The best linear fit to the en¬ 
tire lightcurve is 96.04 counts/d -i- 0.002 counts/d 
HJD, with equal weights for all points. 


We first did a rough period search on the 
delta Ori Aa Chandra lig htcurves using the soft - 
ware package Period04 (ILenz & Bregeri l2005n . 
and found peaks around 4.8d and 2. Id. This 
method uses a speeded-up Deeming algorithm 
(IKurtzl 119851) that is not appropriate for sparse 
datasets such as ours because it assumes the in¬ 
dependence of sine and cosine terms valid only 
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Table 2 

2013 Chandra OBSERVATIONS OF 6 Ori A A 


ObsID 

Start 

HID 

Start 

Phase 

End 

HJD 

End 

Phase 

Midpoint 

HJD 

Midpoint 

Phase 

AT 

days 

Exposure 

s 

Roll 

deg. 

14567 

2456281.21 

396.604 

2456282.58 

396.843 

2456281.90 

396.724 

1.37 

114982 

345.2 

14569 

2456283.76 

397.049 

2456285.18 

397.297 

2456284.47 

397.173 

1.42 

119274 

343.2 

14570 

2456286.06 

397.450 

2456287.52 

397.705 

2456286.79 

397.578 

1.46 

122483 

83.0 

14568 

2456288.67 

397.905 

2456290.12 

398.159 

2456289.39 

398.032 

1.45 

121988 

332.7 


Table 3 

Fourier Periods 


ID^ 

Period 

Amp. 


d 

10“® cts“^ 

raw ctss~^ 

5.0±0.3 

7.1 ±0.7 

residual 

4.76±0.3 

4.7±0.08 

prew. res. 

2.04±0.05 

3.5±0.6 


^raw indicates X-ray lightcurve in ctss“^; 
residual indicates raw lightcurve with linear 
trend removed; prew. res. indicates raw 
lightcuiwe with 4.76d period and linear trend 
removed. 
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Fig. 1.— Chandra X-ray lightcurve from the 2012 campaign with the simultaneous continuous MOST 
optical lightcurve. The time intervals for each of the Chandra observations are delineated with vertical lines 
with the Chandra observation ID (ObsID in Table lljl at the top of the figure. The Chandra lightcurves were 
calculated from the dispersed spectra in each observation. The four observations are separated by gaps due 
to the passage of Chandra through the Earths radiation zone as well as necessary spacecraft thermal control, 
during which time continued 6 Ori Aa observations were not possible. Chandra counts per second are on the 
left y-axis. MOST differential magnitudes are on the right y-axis. 


for regular lightcurves. Therefore, to verify this 
preliminary conclusion and get final resulfs, we 
rafher rely on several mefhods specifically suif- 
able fo such sparse lighfcurves: fhe Fourier pe- 
riod search opfirn a Uy adapfed to spar se datasets 
( Heck et al. 1985 : Cosset et al. 2001, see Fig. 
5), as well as variance and en t ropy methods 
(e.g. Schwarzenberg-Czernvlll989l:ICincotta et al. 


1999h . The results of all these methods were con¬ 
sistent within the errors of each other. 


Using these tools, we first looked for peri¬ 
ods in the raw ctss“^ lightcurve data. A pe¬ 
riod of 5.0±0.3d was found with an amplitude of 
7.1 X 10“^. We then removed the linear trend 
described above from the raw data, producing a 
residual lightcurve. The period searches were re¬ 
peated, with an identified period of 4.76±0.3d 
(amplitude=4.7 x 10“^). After pre-whitening the 


residual data in Fourier space for this period, 
an additional significant period of 2.04±0.05d 
(amplitude=3.5 x 10“^) was found. Each of these 
periods has a Significance Fevel (SF) of SL <1% 
with the definition of SF as a test of the proba¬ 
bility of rejecting the null hypothesis given that 
it is true. If the SL is a very small number, the 
null hypothesis can be rejected because the ob¬ 
served pattern has a low probability of occuning 
by chance. Fig.[3]shows the results of the Fourier 
period search method for the raw, residual, and 
pre-whitened lightcurves, which produced peri¬ 
ods of 5.0d, 4.76d, and 2.04d, respectively. Ta¬ 
ble |3] lists the frequency, period, and amplitude of 
the periods. Fig. |4] shows the residual lightcurve 
with the period 4.76d plotted on the x-axis. The 
residual data were smoothed with a median filter 
for the plot only, in order to see the short-term 
variability more clearly; the analysis used the un- 
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Fig. 2.— Phased Chandra X-ray lightcurve from the 2012 campaign with the simultaneous continuous 
MOST optical lightcurve. The mean 6 Ori Aa MOST optical lightcurve is plotted above the Chandra X-ray 
lightcurve. The four Chandra observations are shown in red, magenta, blue, and black. Binary phase is on 
the x-axis, MOST differential magnitude on the right y-axis, and Chandra count rate on the left y-axis. The 
MOST lightcurves have been smoothed and both lightcurves have been repeated for one binary orbit for 
clarity. 


smoothed residual data points. 

The 5.0d period in the raw data and the 4.76d 
period in the residual data with the linear trend 
removed are considered to be the same period 
because the errors overlap. Comparing the peri¬ 
ods identified in the Chandra 6 Ori Aa data with 
the MOST optical periods, the strongest Chan¬ 
dra period of 4.76d is consistent within the errors 
to the M0ST^2 period of 4.614d (see Table [I])- 
The Chandra period of 2.04d is consistent with 
the less significant MOSTi?io period of 2.138d. 
There is no evidence of an X-ray period matching 
the binary period of 5.73d. 

Finally, we again searched for periods includ¬ 
ing the lightcurve from the original Chandra ob¬ 
servation of 6 Ori Aa (ObsID 639 taken in 2001). 
The lightcurve for obsid 639 alone did not yield 
any statistically significant periods since it cov¬ 
ers a much smaller time interval (60 ks, hence 0.7 


d). However, when combined with the 2013 data, 
we find similar results as mentioned above for the 
raw data, with a period of 5.0d (and its harmon¬ 
ics at 2.5d) providing the strongest peak. Resid¬ 
ual lightcurves were not analyzed with the early 
observation because the linear function (which is 
probably not truly linear) could not be determined 
across an 11-year gap. 

4. TIME-AND BINARY PHASE-RESOLVED 
VARIABILITY 

4.1. Time-sliced spectra 

The discrete photon-counting characteristic of 
the ACIS detector allows the creation of shorter 
time segments of data from longer observations. 
Time-resolved and phase-resolved variability of 
flux and emission line characteristics were inves¬ 
tigated using a set of short-exposure spectra, con- 

















































































Fourier Period Search 



0 0.5 1 1.5 2 

Frequency (day”^) = 1/Period 

Fig. 3.— Periodograms derived using a Fourier period search adapted for sparse datasets. Frequencies 
identified in other wavelengths are shown as red vertical lines. The left red vertical line corresponds to the 
binary period (5.73d), and the two strongest secondary MOST periods (4.613d and 2.448d) are indicated by 
center and right red vertical lines. Top two panels: Periodogram for the raw and residual data, respectively. 
3rd panel: Periodogram for the residual data after “cleaning” (prewhitening) of the strongest signal (4.76d), 
leaving clearly a second period (2.04d). Bottom: Associated spectral window, showing the relative positions 
where aliases may rise. 


for any detailed changes in local response with 
time, such as might be introduced by the Ril ks as¬ 
pect dither pointing of the telescope. This was ac- 


tiguous in time (“time-sliced spectra”), covering 
the entire exposure time of the observations, along 
with individual instrumental responses to account 
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Residual count rate 


0.015 

0.010 

0.005 

0.000 

- 0.005 

-0.010 

0.0 0.5 1.0 1.5 2.0 

phase 

Fig. 4.— The residual count rate lightcurve of 
1 ks binned data after correction for linear fit and 
filtering by a running median over Ril/4d bin size. 
The x-axis is phase for the 4.76d period found us¬ 
ing period search techniques. See text for expla¬ 
nation. 

complished by reprocessing the set of time-sliced 
data using TGCat software, taking care to align 
the zeroth order images among the time-sliced 
spectra prior to the spectral extraction to pro¬ 
duce correct energy assignments for the events. 
The resulting time segments are of similar expo¬ 
sure time, f«40 ks, making them easily compara¬ 
ble. Table|4]lists the 12 time-sliced spectra, along 
with the beginning and end time and binary phase 
range. Our Chandra observations cover parts of 
orbits 396-398 based on the ephemeris. The in¬ 
teger portion of the binary phase is in reference 
to the epoch of the ephemeris used. The decimal 
portion is the binary phase for the specific orbif. 

In addifion fo fhe 12 lime-sliced specfra de¬ 
scribed above, we produced lime-sliced specfra 
of approximately 10 ks in lenglh from fhe Chan¬ 
dra observations, using fhe same fechnique de¬ 
scribed. Forfy-eighl lime-sliced specfra wilh ap¬ 
proximately 10 ks exposure lime each were used 
in fhe composite speclral line analysis in Seel. 

All olher analyses used fhe 40 ks lime- 
sliced specfra. We have no! included fhe 2001 
Chandra observation, obsid 639, in fhe follow¬ 



ing lime-resolved emission-line analyses, primar¬ 
ily because fhe long-lerm frend of flux variabil¬ 
ity seen in our 9-day observing campaign would 
make inferprelafion of fhis early observation ques- 
fionable wilh respeef fo flux level. 

We describe below several differenl analyses 
of fhe variabilify of fhe dispersed speclral dala. 
All comparisons in fhis Section are made fo fhe bi¬ 
nary orbifal period, nof fo fhe periods found in fhe 
X-ray flux in Seef. |3l because we are interested 
in relating any variability fo fhe known physical 
paramelers of fhe syslem and possible effecls of 
fhe secondary on fhe emission from fhe primary 
wind. Firsf, slalisfical fesfs were performed on 
narrow wavelengfh-binned dafa for fhe 12 lime- 
sliced specfra of 40 ks each fo lesl for variabilify 
using fesfs. We Ihen formally til fhe emission 
lines in each of Ihese individual 40 ks lime-sliced 
specfra using Gaussian profiles (Seef. 14.31) . deter¬ 
mining fluxes, line widlhs, and la confidence lim- 
ils. Subsequenl composite line speclral analysis 
used fhe combined H-like ion profiles, as well as 
Fe XVII, fo evaluate fhe flux, velocily, and line 
widlh as a function of binary phase (Secf. l4~4l) . Fi¬ 
nally, we looked for variability in fhe /zr-inferred 
radius (Rfir) of each ion, as well as X-ray lem- 
peralures derived from fhe H-like fo He-like line 
ratios fSecf. l431) . 

4.2. Narrow-band fluxes and variability 

For fhe following slalisfical analysis, we used 
fhe narrow wavelengfh-binned bands in fhe 12 
lime-sliced, 40 ks specfra described above. The 
counl rales for a slandard sel of narrow wave- 
lengfh bins were oulpuf from TGCal processing. 
The paramelers of fhe bins befween 2.5 A and 
22.20 A are lisled in Table |5] We searched for 
variations using a series of fesfs, frying sev¬ 
eral hypofheses (consfancy, linear frend, parabolic 
frend) and checked fhe improvemenf when more 
free pai'amefers were used. The Significance 
Level (SL) is defined in Secl.O Five bands were 
significanfly variable when compared fo a con- 
sfanl value, i.e. SL <1%: (1) fhe continuum 
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Table 4 

Chandra Time-Sliced Spectra Log 


Obsid/slice 

Start HJD 

Stait phase 

End HJD 

End phase Duration (s) 

Mid phase 

14567/1 

56280.718 

396.606 

56281.156 

396.682 

37811 

396.644 

14567/2 

56281.156 

396.682 

56281.607 

396.761 

39000 

396.722 

14567/3 

56281.607 

396.761 

56282.067 

396.841 

39693 

396.801 

14569/1 

56283.267 

397.051 

56283.729 

397.131 

39948 

397.091 

14569/2 

56283.729 

397.131 

56284.204 

397.214 

41000 

397.173 

14569/3 

56284.204 

397.214 

56284.666 

397.295 

39906 

397.254 

14570/1 

56285.568 

397.452 

56286.038 

397.534 

40584 

397.493 

14570/2 

56286.038 

397.534 

56286.524 

397.619 

42000 

397.576 

14570/3 

56286.524 

397.619 

56287.004 

397.703 

41521 

397.661 

14568/1 

56288.177 

397.907 

56288.648 

397.989 

40662 

397.948 

14568/2 

56288.648 

397.989 

56289.134 

398.074 

42000 

398.032 

14568/3 

56289.134 

398.074 

56289.608 

398.157 

40941 

398.116 


Table 5 

TGCat Wavelength Bins 


Label 

A 

A 

A low 

A 

A high 

A 

02500“" 

2.50 

2.00 

3.00 

S XVI 

4.75 

4.70 

4.80 

c4900 

4.90 

4.80 

5.00 

SXV 

5.08 

5.00 

5.15 

c5700 

5.70 

5.40 

6.00 

Si XIV 

6.17 

6.10 

6.25 

c6425 

6.42 

6.30 

6.55 

Si III 

6.70 

6.60 

6.80 

c7800 

7.80 

7.40 

8.20 

Mg XII 

8.40 

8.35 

8.45 

c8800 

8.80 

8.50 

9.10 

Mg XI 

9.25 

9.10 

9.40 

FeXX 

11.20 

10.40 

12.00 

NeX 

12.10 

12.10 

12.20 

cl3200 

13.20 

13.00 

13.40 

Ne IX 

13.60 

13.40 

13.80 

Fe XVII 

15.00 

14.95 

15.05 

cl4925 

14.92 

14.90 

15.05 

O VIII 

16.00 

15.95 

16.05 

C16450 

16.45 

16.20 

16.70 

Fe XVII 

17.07 

17.00 

17.15 

OXIII 

18.98 

18.90 

19.05 

C20200 

20.20 

19.20 

21.20 

0 VII 

21.85 

21.50 

22.20 


‘^Continuum band labels are 
a “c” followed by the band 
wavelength in mA. 
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centered at 4.9 A, (2) S XV, (3) Si XIII, (4) Fe XX 
(I0.4-I2 A), and (5) Ne IX. A further 4 bands are 
marginally variable, i.e. 1% < SL < 5%: (1) 
Mg XII, (2) the continuum centered at 8.8 A, (3) 
Ne X, and (4) the continuum centered at 14.925 
A. When compared to a linear trend, all but Fe XX 
were significantly variable, and when compared 
to a quadratic trend, all but Ne IX, Fe XX, and 
Si XIII were significantly variable, although in all 
cases SL < 5%. 

As an additional test, we directly compared the 
spectra one to another. Using a test on the 
strongest wavelength bins, with spectra binned at 
0.02 A, variability was significant for lines (or in 
the regions of lines): Si XIII, Mg XII, Mg XI, 
Ne IX, and the zone from 10.4-12 A (correspond¬ 
ing to Fe XX). 
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In summary, S XV, Si XIII, Ne IX, and Fe XX 
were variable in both of the above tests. An ex¬ 
ample of the Si XIII lines for several time slices 
is shown in Fig. [51 demonstrating the observed 
variability. As noted later, Ne IX is contaminated 
by Fe lines, which may contribute to the variabil¬ 
ity. A few other emission lines as well as some 
continuum bandpasses may also be variable, but 
with lower confidence. As an additional confirma¬ 
tion of variability for one feature. Si XIII, we per- 
for med a two-sample Kolmogorov Smirnov (KS) 
test ([Press et al.ll2002h on each time slice against 
the complementary dataset (to ensure the datasets 
are independent). With the criteria that an emis¬ 
sion line is variable if the null hypothesis is < 0.1, 
and that a line is not variable if the null hypothe¬ 
sis is > 0.9, the only spectrum where the K-S test 
suggested real variability (at about 2% probability 
of being from the same distribution) was Si XIII 
at (/)=0.11. Note that the KS test shows that there 
is variability, but not what is varying, such as flux 
or cenfroid. 


4.3. Fitting of Emission Lines 

For each of the 12 time-sliced spectra of 40 
ks each, the H- and He-like lines of S, Si, Mg, 
Ne, and O, as well as Fe XVII 15.014 A, were fit 
using the Interactive Spectral Interp retation Sys¬ 
tem (ISIS; Houck & Denicola 20001) . Only Gaus¬ 
sian profiles were considered because (1) Gaus¬ 
sian profiles are generally appropriate at this res¬ 
olution for thin winds at the signal-to-noise level 
of the time-sliced spectra (with some exceptions, 
noted below), and (2) previous studies indicated 
Gaussian profiles provided adequate fits to the 
emission lines for both the early Cha ndra obser¬ 
vation of 6 Ori Aa ( Miller et al. 2002 ) and for the 
combined spectrum from 2012 (Paper I). We note 
that lines in the spectrum of the combined HETGS 
data showed some deviations from a Gaussian 
profile (Paper I). 

The continuum for each time-sliced spectrum 
was fit by using the same 3-temperature APEC 
model derived in Paper I. This model allowed for 


line-broadening and a Doppler shift. Some abun¬ 
dances were fit in order to minimize the residu¬ 
als in strong lines. An Nh of 0.15 x 10^^ cm~^ 
(Paper I) was fixed for the value of the total ab¬ 
sorption. Only the continuum component of this 
model was used for continuum modeling in the 
following analysis. 

Pits were determined by folding the Gaussian 
line profiles through the instrumental broaden¬ 
ing using the RMP and ARP response functions, 
which were calculated individually for each time- 
sliced spectrum. All first order MEG and HEG 
counts, on both the plus and minus arms of the 
dispersed spectra, were fit simultaneously. Por 
most H-like ions, the line centroid, width, and flux 
were allowed to vary. In a few cases where the 
S/N was low, the line center and/or the width was 
fixed to obtain a reasonably reliable fit based on 
the Cash statistic. Por the He-like line triplets, the 
component lines were fit simultaneously with the 
line centroid of the recombination (r) line allowed 
to vary, and the centroids for the intercombina¬ 
tion (/) and forbidden (/) lines forced to be offset 
from the r line centroid by the theoretically pre¬ 
dicted wavelength difference. The individual flux 
and width values of the triplet components for the 
He-like ions were allowed to vary except for a few 
cases of low S/N when the width value of the i 
line and / line were forced to match that of the r 
line to obtain a reasonable statistical fit. The re¬ 
duced Cash statistic using subplex minimization 
was used to evaluate each fit. Por most emission 
lines with good signal-to-noise, the reduced Cash 
statistic was 0.95-1.05. A few of the lines with 
poor signal-to-noise had a reduced Cash statistic 
as low as 0.4. Confidence limits were calculated 
at the 68% (i.e. ztla) level for each parameter of 
each line, presuming that parameter was not fixed 
in the fit. 

In most cases, the line profiles in the time- 
sliced spectra were well fit with a simple Gaus¬ 
sian. In a few cases, a profile might be better de¬ 
scribed as flat-topped with steep wings. Broad, 
flat-topped lines are expected to occur when the 
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Fig. 5.— Si XIII profile (blue) overplotted with the mean Si XIII profile (red) of all time-sliced spectra. 
Upper panel left: Phase is centered at 0.091 and is 0.15 wide; Upper panel right: Phase is centered at 0.254 
and is 0.15 wide; Lower panel left: Phase is centered at 0.576 and is 0.15 wide; Lower panel right: Phase is 
centered at 0.644 and is 0.15 wide. 


formation region is located relatively far from the 
stellar photosphere, where the terminal velocity 
has almost been reached. In such cases, Gaussian 
profiles are expected to fit rather poorly. Occa¬ 
sionally a second Gaussian profile was included 
for a line if a credible fit required it. If more than 
one Gaussian profile was used to fit a line, the to¬ 


tal flux recorded in the data tables is the sum of 
the individual fluxes of the Gaussian components 
with the errors propagated in quadrature. 

For the case of Ne IX where lines from Fe XVII 
and Fe XIX provided significant contribution in 
the wavelength region of the fit, these additional 
Fe lines were fit separately from the Ne IX com- 
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Table 6 

Emission Line Fluxes®' :S and Si. 


Binary Phase 

SXVI 


sxv 


Si XIV 


Si XIII 




r 

i 

f 


r 

i 

f 

0.606-0.682 

1 4+^-^ 
^•^-0.8 

2.6^1;® 

1 9+^-’^ 
^•^-0.9 

1 4+1-8 
^•^-0.8 

2 6+0-9 

^•0-0.7 

ii-e+J 

^•^-0.8 

o.8l?;? 

0.682-0.761 

Q -7 + 0.5 

2.6ti;t 

1 1+^-^ 
^•^-0.8 

1-^-0.8 

3 7+1-1 
■^•'-0.9 

11.31?'? 

4.i1?;®3 

4.5l?;? 

0.761-0.841 

Q -7 + 0.5 

0-3-0.3 

0 9+^0 
^•^-0.6 

0 S+O-”^ 

0 q+2.2 

io.4i?:i 

7-2+:8 

5.0I?:? 

0.051-0.131 

0 1+^-^ 
^•^-0.1 

2 

^•^-1.3 

n 7+1-2 
^+-0.6 

0 7+^0 
^•^-0.6 

5.ill:o 

12.51?'? 

0 + 1-3 
J.J_1 2 

5 4+1-2 
^•^-1.1 

0.131-0.214 

^•°-0.6 

4 2+1-’^ 

*^•^-0.6 

e 2+1-9 
^•"^-1.5 

2.8+:° 

14.2+:? 

^•^-0.4 

0 7+2.0 
^^-2.4 

0.214-0.295 

^•^-0.7 

2 1 +1.5 
J.l_l 2 

1-^-0.7 

^•0-0.9 

4 0 + 1-1 
^•■^-0.9 

14.6+:? 

^•^-0.9 

2 q+0-6 

^•^-0.5 

0.452-0.534 

2C+1.4 

4 2+0-8 

**■3-0.3 

3 . 0 + 5 ^ 

4 6+^’^ 
^•0-1.0 

18.61?:? 

l.O-l.O 

7 2 + 1-8 
'•^-1.2 

0.534-0.619 

^•^-0.6 

4 1+1.8 

^•^-1.5 

2 4+1-^ 
^•^-1.4 

"^•1-1.1 

2 7+1-0 
"^•'-0.8 

Q 7+1-4 
^•'-1.4 

5.01?:? 

6.3l?:? 

0.619-0.703 

^•°-0.6 

4 8+2-1 

T -3 + 1.9 
^••^-1.6 

1.6+;? 

1 2+0.6 

*■■^-0.5 

12.5+:? 

3.7l?:? 

7.4i?:? 

0.907-0.989 

^•^-0.6 

2 7+^-'^ 
^•^-0.5 

**■^-0.8 

2.4l{'t 

2 2+1-2 
^■"^-0.9 

14.6+:? 

^•^-0.9 

5.8l?;? 

0.989-0.074 

^•^-0.7 

2 1 +1-8 
^•1-1.2 

2 4+1-8 
^+-1.5 

^+-0.8 

6.8+:= 

16.01?:? 

2 4+1-3 
^•^-1.2 

9 1 +1-0 

^•^-1.5 

0.074-0.157 

1 

^•°-0.8 

6.8t?;? 

**■*>-0.6 

c q+1-8 

^•^-1.5 

4.4+:? 

19.01?:? 

/•U_1.4 

1.4-1.4 


^10 ^photons!s!cm?. Listed in time order 


Table 7 

Emission Line Fluxes® : Mo and Ne. 


Binary Phase 

Mg XII 


Mg XI 


NeX 


NeIX 




r 

i 

f 


r 

i 

f 

0.606-0.682 

7.4+:? 

25.41?:? 

15.5+? 

9 1 +0-9 

^•1-1.9 

78.3l?i®5 

156.2++ 

77.0I??:? 

16 1+^-9 

0.682-0.761 

6.4++ 

26.1+:? 

14.6l?'? 

9.7i?:? 

66.81?:? 

167.3+?'i 

110.4+?:? 

22.8+'? 

0.761-0.841 

10.31?:® 

22 1 +3.5 
^■^■1-3.2 

1 5 7+8-1 

8.3l?:? 

69.6l?+ 

115.1++ 

152.3+?:? 

28.6+'? 

0.051-0.131 

10.51?:® 

27.0+:? 

18.71?'? 

7.2+:? 

102.81??:? 

145.6+?:? 

147.71??'® 

26 4+7-5 

20.4_g_Q 

0.131-0.214 

lo.ol?:? 

30.41?:? 

8.4+:? 

10.51?:? 

85.3+3=1 

136.3+®:? 

58.31?®/ 

17.61®:® 

0.214-0.295 


28.5+:® 

13.41?:? 

11.2+:® 

81.4l?+ 

196.0l??'? 

64.ll??:? 

1 5 4+7-2 

0.452-0.534 

9.6l?:i 

27.8+:? 

11.8+:? 

10.41®:? 

97 9+19 -1 

213.91??:? 

84 2+18-9 

28.0+'? 

0.534-0.619 

8.8l?:? 

25.71?:? 

1 c 2+2.6 

19-/_2.4 

6.6+:? 

77 

/ /-D-g.S 

147.01??:® 

119.71??:® 

16.o1?®6® 

0.619-0.703 

9.ol?:? 

37.6+:? 

3.6l?:? 

11.61?:? 

70.81?'® 

147 7+19-0 

92.21??:? 

36.9+'? 

0.907-0.989 

10.81?:? 

23.41?:? 

16.51?:? 

11.4+:? 

80.7l?°7® 

177.51??:? 

148.31??:? 

30.1+'? 

0.989-0.074 

9 2+^-'^ 

22 2 + 2.3 
^"^■■^-2.9 

14.91?'® 

8 9+^-^ 

92.9l?+ 

175.5+?:® 

138.41?®:® 

43.7+:? 

0.074-0.157 

12 3+^-® 

30.21?:? 

15.5+:? 

9 4+2-8 
^•^-2.0 

69.2l?°8= 

157.5+®:? 

61 . 21 ??:? 

22 2 + 7-2 

zo.z,_6 2 


® 10 ^photons/s/cm?. Listed in time order 
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Table 8 

Emission Line Fluxes®' : O and Fe XVIL 


Binary Phase 

0 VIII 


0 Vll 


Fe XVII 



r 

i 

f 


0.606-0.682 

1018.6/®3.i 

1101 

841.2/1®®;® 

QO q75.2 
^^•^51.9 

299.0/=/l 

0.682-0.761 

QQ9 1 +91-9 

408 

1058.2/=®/® 

q65.4 

92.6347 

304.3/1/,? 

0.761-0.841 

1010 

Q1 T £+208.5 
j.u_2Q2.3 

889.8/“/= 


270.7/1/® 

0.051-0.131 

1023.7l92| 

£££ 1+152.1 

666 . 1 _j 3 j g 

/VD.6_i47 2 

13.3i5,3 

307.8/=/,! 

0.131-0.214 

1122.8/9^;® 

84:5.9_igg 0 

/U‘t.i_i5i 0 


243.i/=/®5 

0.214-0.295 

920.4/«^;S 

728.1/1^// 

650 

D9U.D_j^3g g 

13.3i5,3 

275.8/=// 

0.452-0.534 

1018.6/^;® 

990.0/=®// 

79T c + i97.7 

/Z9.J_j^gg g 

13.435 4 

310.2/=/=5 

0.534-0.619 

1066.0l«/ 

880.9/1®/,® 

586.3/«/® 

21.921.9 

294.5/=/,® 

0.619-0.703 

860.9/®®/ 

1215.4/1®/,® 

398.6/1=/® 

25.8®®:1 

254.7/=/® 

0.907-0.989 

1022.4/®/® 

689.3/14/® 

738.6/1®/® 

lI-2f/J 

9U9.2_25 6 

0.989-0.074 

1059.4/®/® 

683.7/®!®-® 

750.0/1®/® 

26.5|/g 

321 . 0 /=/® 

0.074-0.157 

1026.2/®/® 

844 5+157.8 

644.9 _j4Q 7 

684.8/11/® 

6.6|®6® 

375.0/=®;! 


® 10 ^photons/s/cm?. Listed in time order. 


ponent lines. The Fe lines fit in this region were 
Fe XIX at 13.518 A, Fe XIX at 13.497 A, and 
Fe XVII at 13.391 A, with Fe XIX at 13.507 A 


and Fe XIX at 13.462 A included at their theoret¬ 
ical intensity ratios to Fe XIX at 13.518 A. Also, 
the Ne X line is blended with an Fe XVII line. We 
have assumed this Fe XVII component contributes 
flux to the N e X line equal to 13% of the Fe XVII 
at 15.014 A ( Walborn et al.ll2009h . A final correc¬ 
tion was applied to the Si XIII-/ line because the 
Mg Lyman series converges in this wavelength re¬ 
gion. Using theoretical relative line strengths, we 
assumed the Si XIII-/ line flux was overestimated 
by 10% of the measured flux of the Mg XII La 
line. 


The flux values and confidence limits are tabu¬ 
lated in Table 14.31 for S and Si lines. Table 14.31 for 
Mg and Ne lines, and Table [8] for O and Fe XVII 
lines. To summarize the results. Fig. [6] shows a 
comparison of the fluxes of the H-like ions. The 
error bars for S XVI are quite large. Si XIV shows 
a peak at about />=0.0 and a somewhat lower value 
at about (j)=0.6. Mg XII, Ne X, and O VIII are 
essentially constant. Fig. |7] shows the fluxes for 


the He-like r lines. S XV-r has a maximum at 
about (/)=0.1, with lower flux in the range i/=0.5- 
0.8. The flux values for Si Xlll-r are larger for 
/=0.0-0.4 than for the range /=0.5-0.8. O Vll-r 
shows an apparent increase in flux in the i/=0.5- 
0.7 range. Mg Xl-r and Ne IX-r are relatively 
constant. Ne IX was consistently variable in the 
narrow-band statistical tests. In this Gaussian fit¬ 
ting of the lines we have fit and removed the con¬ 
taminating Fe lines in Ne IX, possibly removing 
the source of variability in this triplet. We note 
that the points in Fig. [6] and Fig. |7] are from three 
different orbits of the binary. The increase in flux 
with time discussed in Sect. [3] has not been taken 
into account in these fitted line fluxes, either in the 
plots or the data tables, so care must be taken in 
their interpretation. 

4.4. Spectral Template and Composite Line 
Profile Fitting 

In order to improve the signal-to-noise ratio in 
line fits in the time-sliced data, we used two meth¬ 
ods to fit multiple lines simultaneously. In the 
first method, we adopted a multi-thermal APEC 


16 










Flux (photons/s/cm^2) 


4x10"®F 


3x10"® : 


0.5 1.0 


OL 

0.5 


01 _, 

0.5 


1.5 

Phose 

perlostron 


1.5 

Phase 


S XVI 


Mg XII 


1 


^ 2X10“' - 


01 _, 

0.5 


1.2x10 ; 

1.0X10““ - 
8.0x10“' - 
6.0x10“' - 
4.0x10“' - 
2.0x10"' - 


OL 

0.5 


1.5 

Phose 

perlostron 


Si XIV 


Ne X , 


1.0 1.5 2.0 

Phase 


0 VIII 


1.5 

Phase 


Fig. 6.— Fluxes of H-like emission lines based on Gaussian fits vs phase. Errors are la confidence limits. 
Phase with respect to periastron is indicated at the top of the plot. 


Smith et al. 

2001 


Eoster et al. 

2012) 


model which describes the mean spectrum fairly 
well (see Table 14.41) . and used this as a spectral 
template, allowing the fits to the Doppler shift, 
line width, and overall normalization to vary 
freely. This is a simpler model than the more 
physically based APEC model defined in Paper 
I since here it need not tit the spectrum globally, 
but is only required to tit a small region about 
a feature of interest. To demonstrate this, Eig- 


ure[8]shows an example after fitting only the 8.3- 
8.6 A region for the Mg XII feature’s centroid, 
width, and normalization for the entire exposure. 
The temperatures were not varied, and the rela¬ 
tive normalizations of the three components were 
kept fixed, as was the absorption column. We can 
see that this provides a good local characteriza¬ 
tion of the spectrum, and so will be appropriate 
for studying variation of these free parameters in 
local regions as a function of time or phase. Eor 
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Fig. 7.— Fluxes of He-like r emission lines based on Gaussian fits vs phase. Errors are la confidence 
limits. 


any such fit, other regions will not not necessarily 
be well described by this model. 

These fits used spectra extracted in 10 ks time 
bins (about 0.02 in phase), but were fit using a 
running average of three time bins. We primarily 
used the H-like Lyman a lines, as well as some 
other strong and relatively isolated features (see 
Table 14.41) . The results for the interesting parame¬ 
ters, the mean Doppler shift of the lines and their 
widths, are shown in Fig.|9l The smooth curve in 


the top panel is the binary orbital radial velocity 
curve. We clearly see a dip in the centroid near 
(/> = 0.8, as well as significant changes in average 
line width. 

In the second method, we rebinned regions 
around selected features to a common velocity 
scale and summed them into a Composite Line 
Profile (“CLP”). While this mixes resolutions 
(the resolving power is proportional to wave¬ 
length and is different for HEG and MEG) and 
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Table 9 

Plasma model parameters used for spectral template fitting. 


Temperature 

Components^ 

T 

Norm 

2.2 

8.16 

6.6 

1.90 

19.5 

0.226 

Relative Abundances'^ 

Elem. 

A 

Ne 

1.2 

Mg 

0.7 

Si 

1.6 

Fe 

0.9 

Total Absorption^ 

Nh 

0.15 


^Temperatures ai‘e given in 
MK, and the normalization is 
related to the volume emis¬ 
sion measure, VEM, and 
distance, d, via VEM = 
10^^ {And^) X EiNorrui 

'^We give elemental abun¬ 
dances r elative to the solar 
values of lAnders & Grevessd 
for those significantly 
different from 1.0. (These 
are not rigorously determined 
abundances, but related to 
discrete temperatures adopted 
and actual abundances.) 

"^The total absorption is 
given in units of 10^^ cm~^. 
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Table 10 

Lines used in CLP analysis'^ 


'^0 

A 

Feature 

6.182 

SiXIV 

8.421 

Mg XII 

10.239 

NeX 

11.540 

Fe XVIII 

12.132 

NeX 

14.208 

Fe XVIII 

15.014 

Fe XVII 

15.261 

Fe XVII 

16.005 

Fe XVIIU-O VIII 

16.780 

Fe XVII 

17.051 

Fe XVII 

18.968 

O VIII 


^ines used in the ensemble fit¬ 
ting with Composite Line Pro¬ 
file or spectral template meth¬ 
ods. Both HEG and MEG were 
used at wavelengths shorter than 
16 A. The region widths were 
0.20 A, centered on each feature. 


blends in the CLP, the mix should be constant 
with phase and be sensitive to dynamics as long 
as line ratios themselves do not change. Hence, 
we can search for phased variations in the line 
centroid. This technique has been applied fruit- 
fully in characteriz i ng steU ar activity in cool stars 


(IHoogerwerf et ahl l2004l : iHuenemoerder et al. 


2006h . The CLP profiles were computed in phase 


bins of 0.01, but grouped by 5 bins for fitting, 
thus forming a running average. We used the 
same lines as in the template fitting. In Fig.fTOlwe 
show an example of composite line profiles, and 
fits of a Lorentzian (since the composite profile 
is no longer close to Gaussian) plus a polynomial 
to determine centroid and width. This method, 
while less direct than template fitting, did con¬ 
firm the trend seen in line velocity in the template 
fitting. 


The template-fit line width result is very in¬ 
teresting in that it shows a significantly narrower 
profile near 0 fhan at other phases. Given the 
trends in width, (low near phase 0.0, high near 0.2 


and 0.8) 10 ks time-sliced spectra were grouped 
in these states and then compared (Fig. fTTI) . The 
plots show the narrow state in blue and the broad 
state in red. The lines are all sharper, except for 
the line at 17 A (and maybe Si XIV 6 A), in the 
narrow state. The top panel of Fig. [H] shows 
a heavily binned overview, and the lower panel 
shows a compaiison of the Ne X line profiles at 
(/)=0 and at quadrature phases. 

The line width variability was confirmed by 
comparing the average spectrum at phases near 
(j) = 0.0 with the average spectrum at other 
phases. The changes were primaiily in a reduced 
strength of the line core in the phases when the 
lines are broad, with little or no change in the 
wings. 

Fig. [12] shows the trend vs. emission line. Ex¬ 
cept for Ne X and Fe XVll 17 A, there is a 
trend for larger differences in the FWHM with 
increasing wavelength. Note that temperature 
of maximum emissivity goes roughly inversely 
with wavelength; wind continuum opacity in- 
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Fig. 8.— A portion of the HEG spectrum for the 
entire exposure, after fitting an APEC template 
to the 8.3-8.6 A region. The temperatures, rel¬ 
ative normalizations, and absorption were frozen 
parameters. The Doppler shift, line width, and 
normalization were free. We show the resulting 
model evaluated over a bit broader region that was 
fit to demonstrate the applicability of the model to 
the local spectral region. Other regions will not 
necessarily be well represented by the same pa¬ 
rameters. 

creases with wavelength. The increasing trend 
is typical of winds, since the opacity causes 
longer wavelength lines to be weighted more to 
the outer part of the wind where the velocity is 
higher. Gaussian fit centroids show the lines are 
all slightly blueshifted, which could be consistent 
with skewed wind profiles. The “narrow” group is 
near the primary star radial velocity shift. Radial 
velocities of the lines are all roughly consistent 
with -80 km s“^, except perhaps Ne X, which is 
blended with an Ee line. The dependence between 
line width and binary phase was confirmed inde¬ 
pendently by moment analyses of the individual 
lines, and was also suggested by the CEP analysis 
above. 



Pig. 9.— The mean emission line Doppler veloc¬ 
ity (top: points with en'or bars), primary radial 
velocity (top, sinusoidal curve), and mean line 
width for Ne X (bottom) derived by fitting spectra 
in phase bins with an APEC template, allowing 
the Doppler shift, line width, and normalization 
to vary freely. Data from the individual Chandra 
observations are differentiated with colors. Error 
bars (Icj) are correlated over several bins since a 
running average was used over 3-10 ks bins. 


(forbidden, intercombination, resonance) lines by 
defining two rat io jl: the R ratio = f /i an d the G 
ratio=(i + f)/r. iGabriel & JordanI ( 1969lf demon¬ 
strated that the f /i and G ratios are sensitive to 
the X-ray electron density and temperature, re¬ 
spectively. These ratios have been used exten¬ 
sively in stellar X-ray studies. In addition, the 
presence of a strong UV/EUV radiation field can 
change the interpretation of the f /i ratio from a 
density diagnostic to a measurement of the radi¬ 
ation field geometric dilution factor, i.e., effec¬ 
tively the radial location of th e X-ray emission 


from a central radiation field (iBlumenthal et al. 


19721) . The f /i ratio is known to decrease in the 


4.5. X-ray Emission Line Ratios 

The He-like ions provide key plasma diag¬ 
nostics using the relative strengths of their fir 


^ These lines are often designated as w, x, y, z in order to 
highlight the fact that the i-line emission is produced by 
two transitions (x-l-y) such that R = zjG + y) and G = 
{x-\-y-\- z)/w 
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Fig. 10.— An example of Composite Line Pro¬ 
files for two different phases with different cen¬ 
troids, 0.68 (top) and 0.78 (bottom), as defined by 
the photometric ephemeris. In each large panel, 
the histogram is the observed profile, fhe smooth 
curve is the fit. In the small panel below each are 
the residuals. 


case of a high electron density and/or high radia¬ 
tion flux density, which will de-populate the upper 
level of the /-line transition (weakening its emis¬ 
sion) while enhancing the f-line emission. For hot 
star X-ray emission, the /// ratio is controlled en¬ 
tirely by the strong UV/EUV photospheric radia¬ 
tion field. The first analysis of an O supergian t 
HETG spectrum by Waldron & Cassinelli ( 2001 ) 
verified that the observed X-ray emission is dis¬ 
tributed throughout the stellar wind and demon¬ 
strated that density effects could only become im¬ 
portant in high energy He-like ions if their X-rays 
are produced extremely close to the stellar sur¬ 
face. Thus the /// ratio can be exploited to deter¬ 
mine the onset radius or fir-inferred radius {R fir 
in units of R*) of a given ion via the geomet¬ 
ric di lution factor of the photosp heric radiation 
field ( Waldron & Cassinelli 2001 ). In addition. 



Wavelength (A) 



Eig. 11.— Examples of broad and narrow emis¬ 
sion lines for selected wavelength regions. Plots 
are constructed from counts per bin data with¬ 
out continuum removal. Eor this comparison, the 
10 ks time-sliced spectra have been combined to 
represent (/=0.0 (blue) and the quadrature phases 
(red). 


there are basically two types of /ir-infeiTed radii, 
“localized” (point-like) or “distributed.” The first 
detailed distribute d appr oach was developed by 
Eeutenegger et al. ( 20061) assuming an X-ray op¬ 


tically thin wind. Eor a given observed f /i ra¬ 
tio, the localized approach predicts a larger Rfir 
as compar ed to the distributed approac h (see dis¬ 


cussion m 


Waldron & Cassinellill2007l) . Since all 


X-ray emission lines scale as the electron density 
squared, all line emissions are primaiily domi¬ 
nated by their densest region of formation. How¬ 
ever, in the case of the fir lines, an enhanced i- 
line emission can only occur deep within the wind 
(high density), whereas the majority of the /-line 
emission is produced in the outer wind regions at 
lower densities. The r-line emission is produced 
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Fig. 12.— Comparison of FWHM in km s“^ 
for several emission lines in the time-sliced spec¬ 
tra of S Ori Aa. For this comparison, the 10 ks 
time-sliced spectra have been combined to repre¬ 
sent (/)=0.0 (blue) and the quadrature phases (red). 
Note that these spectra do not have continuum re¬ 
moval. Gaussian plus polynomial line fitting was 
used on the time slices to determine the line width. 


throughout the wind. 

Another X-ray temperature-sensitive line ratio 
is the H-like to He-like line ratio ( H/He) as ex 


2002 

: Miller et al. 2 

o 

o 

Waldron et al. 20041: 

Waldron & Cassinelli 

2007 

). However, a wind 
implies a density de- 
nd He-like lines may 
ns) and a dependence 
orption effects. Thus, 
m H/He ratios mav be 

distribution of X-ray sourcef 
pendence (i.e., the H-like a 
be forming in different regie 
on different wind X-ray abs 
the temperatures derived fro 
higher than their actual valu 

es tsee Waldron et al. 


20041 : IWaldron & Cassinellill2007l) . 


Our line ratio analysis is based on t he ap - 
proach given by IWaldron & Cassinellil (120070 . 
The f/i ratios for each He-like ion in each 
time-sliced spectrum are tabulated in Tables fTTl 
-M We did not include S in this analysis be¬ 
cause the flux measurement errors are large and 
the flux ratio errors are extremely large or un¬ 
bounded. We calculated the /zr-inferred radii 


(Rfir) and H/He-inferred temperatures (Thhe) 
versus phase for the 12 time-sliced spectra, as de¬ 
termined by the Gaussian line fitting. All radii 
were determined by the point-like approach and 
a TLUSTY photospheric radiation field with pa¬ 
rameters Te//=29500 kK and log G=3.0. The 
model f/i ratios and H/He ratios used to extract 
Rfir and Thhc information take into account the 
possible contamination from other lines. For all 
derived Rfir, we assume that the He-like ion line 
temperature is at its expected maximum value. 

The /zr-inferred radii for Mg and Si are plot¬ 
ted in Fig. [13] and for Ne and O in Fig. [T4| In 
all cases the derived Rfir is based on the average 
of the minimum and maximum predicted values, 
so if the lower bound is at 1, then upper bound 
should be considered as an upper limit. In these 
plots, any lower limits are indicated by arrows. 
The binary phase is used for the x-axis. There are 
10 cases showing a finite range for the O Rfir- 
There are two O and Ne Rfir values at phase 

0.65, but in different binary orbits, indicating 
the same finite radial locations (within the errors) 
for each ion (O at Ri 7 — 9, and Ne at « 4 — 6), 
which could mean that at least for O and Ne the 
behavior is repeatable. This is not seen in Mg or 
Si. For Mg, in one case for (f> R30.65 there is a fi¬ 
nite Rfir (~2), whereas the other case at (f> RiO.65 
indicates only a lower bound of Ri4.5. For Si 
there are five Rfir with finite ranges, all within 
the errors of one another. Si has four Rfir at 1 
since the observed f/i for these phases were be¬ 
low their respective minimum f/i. This behavior 
suggests that these regions producing the majority 
of the higher energy emission lines may be experi¬ 
encing significant dynamic fluctuations in density 
and/or temperature. 

The 12 derived H/He temperatures {Thhe) 
versus phase for each ion are shown in Fig. [TS] 
In general, for all ions there is very little variation 
in Thhc with phase, although one could easily 
argue that at certain phases there are minor fluc¬ 
tuations. 

A verification of these results was obtained 
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using the Potsdam Wolf- Rayet (PoWR) code 
(IHamann & Grafenen l2004h to perform a similar 
f /i analysis that included diffuse wind emission 
and limb darkening. Differences in the results 
obtained using the two methods were negligible 
compared to measurement uncertainties. 



0.0 0.2 0.4 0.6 0.8 1.0 

phase 


Fig. 13.— Phase dependence of the derived Mg 
and Si /ir-inferred radii (Rjir) for the 12 time- 
sliced spectra. Upper and lower limits are shown 
as arrows. See text for model details. 

4.6. Non-detection of Stellar Wind Occulta- 
tion Effects 

One goal of this program was to use the vari¬ 
able occultation of the primary wind by the essen¬ 
tially X-ray-dark secondary, 5 Ori Aa2, as it or¬ 
bits the primary, mapping the ionization, temper¬ 
ature, and velocity regimes within the primary’s 
stellar wind. The secondary star, 6 Ori Aa2, is or¬ 
biting deep within the wind of the primary star, 
6 Ori Aal. Based on the binary separation of 2.6 



0.0 0.2 0.4 0.6 0.8 1.0 

phase 


Fig. 14.— Phase dependence of the derived O and 
Ne /ir-inferred radii (Rfir) for the 12 time-sliced 
spectra. Upper and lower limits are shown as ar¬ 
rows. See text for model details. 


R^ai and our calculations of the /ir-inferred radii 
of the various ions, we expect the secondary to 
be outside the onset radius of S emission in the 
primary wind, very close to or inside of the on¬ 
set radius of Si, and inside the onset radii of Mg, 
Ne, and O. We can make a simple model of the 
expected lightcurves for these emission lines for 
6 Ori Aa. If the secondary is outside the onset ra¬ 
dius of an ion, the lightcurve will have a maxi¬ 
mum and relatively constant flux value between 
4> ~0.25 and ?«0.75. The lightcurve will have a 
relatively constant but lower flux at / R:^0.75-0.25 
as it occults both the back and front sides of the 
onset-radius shell. If the secondary is inside the 
onset radius of an ion, then the lightcurve will be 
at a maximum flux near 0= 0.0 and 0.5. Between 
these two phases, the lightcurve will have a lower 
and relatively constant flux value as it occults only 
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Fig. 15.— Phase dependence of the Si, Mg, Ne, 
and O Thhc calculated from the H/He ratios of 
each of the twelve 40 ks time-sliced spectra. 

the back side of the onset-radius shell. 

The i line of the He-like triplet is formed deep 
in the wind while both the r and / lines are more 
distributed throughout the wind. Our best chance 
of identifying occultation effects is probably from 
the fluxes of the i lines, as the shell of i-line emis¬ 
sion will be thinner than either the /-line emis¬ 
sion shell or the r-line emission shell. To estimate 
any occultation effect that we might see in these 
lightcurves, we calculated the maximum volume 
of z-line emission that could be occulted by the 
secondary using various parameters. The total 
i-line emission expected from an emitting shell 
around the surface of the star was estimated using 
the spherical volume of the shell. The secondary 
star is assumed to have a radius of 0.3 Rstar- The 
percentage of i-line emission occulted by the sec¬ 
ondary star will be maximum when the largest 


column of emitting material is occulted, which is 
at the rim of the shell. An estimate of the occulted 
emission was based on the volume of the spheri¬ 
cal cap of the emitting shell, with a height equal 
to the diameter of the secondary star. Care was 
taken to subtract the volume of the star that might 
be included in the cap. We find that, for the ex¬ 
treme case of an i-line emission shell at the sur¬ 
face of the primary star with a thickness of 0.01 
R^ai. the maximum occultation of the i-line flux 
would reduce the flux by w20%. If the thickness 
of the Aline shell is instead a more likely value 
of 0.1 R^ai but still in contact with the surface 
of the star, the flux of the i emission line would 
be reduced by only about 3% due to occultation 
by the secondary. In these estimates, the sec¬ 
ondary would not be at primary minimum (^=0.0, 
but instead projected on the rim of the Aline emis¬ 
sion shell where the column of Aline emission is 
greatest. A thicker i-line shell would reduce the 
amount of occultation further because of the finite 
size of the projected secondary star in relation to 
the volume of the sphere of Aline emission. Any 
other position of the secondary in the orbit, or any 
larger shell of Aline emission, would reduce the 
percentage of occultation. 

Fig. [T^ shows the flux measurements for the i 
lines of the He-like triplets. We have previously 
found a linear increase in flux over the time pe¬ 
riod of the observations which is not removed in 
these plots. While the estimates of the fluxes from 
the S XV-i and Si Xlll-i lines that might have an 
onset radius very close to the stellar surface of the 
primary star do not preclude the existence of oc¬ 
cultation effects, we unfortunately cannot identify 
such variability which would be at the 1-2% level, 
particularly due to the other identified variations 
on the order of 10-15% and the errors of the mea¬ 
surements. 
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Fig. 16.— Flux of He-like i component of the 
triplet based on Gaussian fits vs phase. Errors are 
Icr confidence limifs. Fluxes have been normal¬ 
ized by dividing by fhe mean flux for fhe respec¬ 
tive line. Phase wifh respecf fo periasfron is indi¬ 
cated af fhe fop of fhe plof. The plof for each ion 
is offsef by a value of 3 from fhe previous one for 
clarify. 

5. DISCUSSION 

5.1. Effects of a Wind-Star Collision 

An important and unexpected result of this 
analysis is the discovery of variability of line 
widths with binary phase. H-like emission line 
widths are at a minimum at (^=0.0 when the sec¬ 
ondary is in front of the primary, and to a lesser 
degree at (/)=0.7. The line widths are at a maxi¬ 
mum near 0= 0.2 and (/)=0.S, close to quadrature. 
The phase-dependent variability of the emission 
line widths must therefore be related to interac¬ 
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tion between the primary and the secondary. 

In Paper 1 we developed a model to represent 
the effect of the secondary star on the wind region 
of the primary. Our ID, line-of-centers, CAK cal¬ 
culations presented in Paper 1 showed that radia¬ 
tive braking does not occur in this system, so the 
primary wind directly impacts the surface of the 
secondary star. A similar example has been ob¬ 
served for CPD —41°7742, along with an eclipse 
of the X-ray emitting colliding wind r egion when 


the tw o stars were perfectly aligned (IS ana et al 


20051). A 3 D smoothed pa r ticle hydrodynamics 


(SPH) code (iRusselll l2013l : iMadura et al.l 1201 3h 


was then used to simulate the effect of the wind- 
wind collision in 6 Ori Aa. The colliding winds 
form an approximately cone-shaped cavity in the 
wind of the primary with the secondary star at the 
apex of the cone. The cavity has a half-opening 
angle of Ri30°, so the solid angle fraction taken up 
by the cone is «8%. A bow shock suiTounds the 
cavity, and within this cavity the secondary wind 
prevails, yielding lower densities as well as very 
little X-ray flux. Fig.[l7l reproduced from Paper 1, 
shows the density and temperature structure of the 
winds and their interaction in the binary orbital 
plane. The hot gas in the shock at the interface 
between the lower density cavity and the primary 
star wind was calculated in Paper 1 to produce at 
most 10% of the observed X-ray flux. Accord¬ 
ing to this model, as the cavity rotates around the 
primary star from the blue- to red-shifted part of 
the wind, emission line profiles should change in 
shape and velocity (as long as the radius of line 
formation is similar to or larger than the location 
of the apex of the bow shock). 


The variability of the emission-line widths we 
have observed may potentially be explained by 
this cavity in the primary wind caused by the wind 
interaction with the secondary. When viewed at 
(f)=0.0, the cavity will occupy a region of the pri¬ 
mary stellar wind that would otherwise be the 
formation region of emission with high negative 
velocities. The emission line profiles viewed at 
this phase might then be truncated at the largest 
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Fig. 17.— Model of density and temperature structure of the binary orbital plane of the SPH simulation 
of Aal (larger black circle) and Aa2 (smaller black circle), based the parameters from Paper IV and the 
model in Paper I. The collision of the wind of the primary star against the secondary star produces a low 
density cavity within the primary wind. The perimeter of the low density cavity is a shocked bow shock 
of higher density than either star’s wind region (left panel). In the temperature plot on the right, only the 
hot gas from along the wind-collision, bow shock boundary is shown since the SPH simulation does not 
include the embedded wind shocks of either delta Ori Aal or delta Ori Aa2. The primary wind collides 
directly with the secondary surface in this simulation primarily because of the large difference in mass loss 
rate {MAai/MAa 2 ~ 40). 


negative velocities, creating a comparatively nar¬ 
rower profile than one expects without the pres¬ 
ence of the secondary. At (/)=0.5, some of the 
most positive velocities would not be detected in 
the emission-line profiles. Wind absorption must 
be taken into account at (/)=0.5 and the effect of 
the cavity may be less pronounced. Such emis¬ 
sion truncation at (/)=0.0 and 0.5 is suggested in 
Fig. [TSl displaying the Mg XII profiles and fits 
for the time-sliced spectrum near (/)=0.0 and the 
time-sliced spectrum near (/)=0.5. The Mg XII line 
was chosen as a relatively strong line that was also 
used in the template analysis in Sect. 14.41 The red 
line in the figure is the Gaussian fit to the line, 
the black line is the time-sliced spectrum for a 
particular phase interval, and the lower panel of 
each plot shows the 6x statistic (though the Cash 
statistic was actually used in the fitting). Negative 
velocities appear somewhat under-represented in 
the time-sliced spectrum near (/)=0.0, although the 


profile near i;^=0.5 does not appear to be asym¬ 
metric. There could be other explanations for 
this characteristic, such as velocity changes in 
the centroid or non-Gaussian profiles. At quadra¬ 
ture, (/)=0.25 and 0.75, the velocities normally 
produced in the embedded wind of the primary 
and replaced by the region occupied by the cav¬ 
ity will tend to be near zero velocity, resulting 
in emission lines of expected widths, but some¬ 
what non-Gaussian peaks, such as flat-topped or 
skewed peaks. This prediction is not inconsistent 
with the profiles observed near quadrature in the 
6 Ori Aa time-sliced spectra, which are broader 
than the profiles seen near conjunction (see Fig. 
[18] as an example). 

However, furfher analysis of fhe Mg XII line 
reveals a more complex scenario. Fig. [T9| shows 
the correlation of FWHM with flux for Mg XII for 
the 12 time-sliced spectra, based on the Gaussian 
fits. While most of the points show a similar dis- 
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Fig. 18.— Mg Xll profile overplotted with the Gaussian fit. Upper panel left: Phase is centered at 0.031; 
Upper panel right: Phase is centered at 0.492; Lower panel left: Phase is centered at 0.254; Lower panel 
right: Phase is centered at 0.721. 


tribution with little or no dependence of FWHM 
on flux, the three points from the first Chandra ob¬ 
servation, ObsID 14567, have significantly larger 
FWHM values along with lower flux. We have 
carefully checked that there is no known reason 
to expect this result to be instrumental. We con¬ 
clude that the variability in the Mg Xll line has 


several time scales, of which the binary orbit and 
colliding winds associated with the binary system 
are only one component. 
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Fig. 19.— Gaussian parameters for Mg XII lines 
for each of the 12 time-sliced spectra. Error bars 
are Icj. Flux vs FWHM. 


5.2. Stellar pulsations and Corotation Inter¬ 
acting Regions 

Pulsations have been found in only a few O 
stars, but those that have been found are mostly in 


late-0 stars, such as ( Oph, 09.5V (IWalker et al 


even 

over long timescales (i.e. Prin 

W - -, 

a & Howarth 

1986; 

Kaoer et al. 

1996. 1997.1999 

: iMassa et al. 

1995; 

Prinia 1988 

: Henrichs et al. 19881. Models 


predict that massive stars near the main sequence 
will experience pulsations of several types, in¬ 
cluding non-radial pulsations (N RP), 0 Cep insta¬ 


bilities, and Z-mod e pulsations (ICox et al.lll992l: 
Pamyatnykh 1999h . 


Another source of variability relates to phe¬ 
nomena near or on the surface of the star, such 
as bright spots, which tend to be transient over a 
few stellar rotations or less. For many years DAC 
vaiiability in UV P-Cygni profiles has been rec- 
ognized, and is believed to be common in O stars. 
Cranmer & Owocki ( 1996h provided a model of 


ad hoc photospheric perturbations in the form of 
bright spots that was able to reproduce the DAC 
phenomenon. CIRs, possibly related to DACs, are 
perturbations that start at the base of the stellar 


wind and can extend far out into the wind, and 
thus are tied to the rotation period of the star, al¬ 
though they could also be transient. Presumably 
there can be multiple CIRs distributed over the 
surface of the star, cumulatively resulting in vari¬ 
ations that appear to be shorter than the rotation 
period, although each individual spot and its as¬ 
sociated CIR will rotate with the star. 

The periods we have identified in the Chan¬ 
dra lightcurve of 6 Ori Aa (4.76d and 2.04d, see 
SectU, are possibly the X-ray signature of CIRs 
or pulsations. The 4.76±0.3d period may be the 
rotation period of the primary star. Based on a 
V sin i value of 130 km s“^ (Paper IV) and the es¬ 
timated value of R*= 12 R 0 , and also including the 
assumption of alignment between the rotational 
and orbital axes, the rotational period should be 
about 4.7d, consistent with our strongest period 
of 4.76d. A single non-transient CIR would share 
the stellar rotational period. We did not find evi¬ 
dence of the binary orbital period of 5.73d in the 
X-ray lightcurve, indicating that most of the X- 
ray variability is not related to the orbital motion, 
at least for the limited orbital data we have. 

There is an apparent increase in flux over 9 
days, suggesting long-term variability. We do not 
have a sufficient time baseline to quantify this 
component of variability, but it is clearly not as¬ 
sociated with the binary period. We suggest this 
long-term variability is due to pulsations such as 
NRP and/or increasing CIR activity. 

6. CONCLUSIONS 

Chandra high-resolution grating X-ray spectra 
of (5 Ori Aa acquired in 2012 with a total exposure 
time of Ri479 ks have been analyzed for phase- 
resolved and time-resolved variability. Several 
components of variability were detected in our 
analyses. 

The count rate of the entire spectral range in¬ 
creased during the 9-day observing campaign by 
approximately 25%. We cannot constrain the 
cause of this longer-term vaiiability with this 


29 






































































dataset, but we speculate this may be related to 
stellar pulsations or CIRs and other wind instabil¬ 
ities. 

An important result of the period searches in 
the X-ray data is the binary motion seems to con¬ 
tribute very little to the variability of the total flux. 
A period search of the total X-ray flux lightcurve 
yielded periods of 4.76±0.3d and 2.04zb0.5d af¬ 
ter removal of the long-term trend, both of which 
are less than the binary period of 5.73d. A pe¬ 
riod search including the early 2001 Chandra ob¬ 
servation as well as the 2012 data gave a period 
near 5.0d, within the errors of the 4.76d period de¬ 
termined from the normalized 2012 spectra. The 
4.76d period is consistent with the secondary pe¬ 
riod found by MOST of 4.614d; thus it is present 
in both X-ray and optical data. We suggest that 
this may be the rotation period of 6 Ori Aal based 
on estimates of v sin i and the radius of 6 Ori 
Aal. The 2.04d period, also found in the MOST 
photometry, may be associated with pulsations or 
CIRs. 

Flux variability of individual emission lines 
was confirmed with statistical tests for the He- 
like triplets of S XV, Si XIII, and Ne IX (con¬ 
taminated with an Fe XVII line), as well as the 
Fe XX complex. Also, several line profiles are ap¬ 
parently non-Gaussian with blue-shifted centroids 
of about -80 km s“^ prevalent, possibly indicating 
that line-fitting with wind profiles would be more 
appropriate. Derived Rfir are in similar ranges 
for O stars of the spectral type of 6 Ori Aal. 

For first time, phase-dependent variability in 
the X-ray emission line widths has been found in 
a binary system. Line widths are at a minimum 
at 0=0.0 and at a maximum at 0=0.2 and 0.8, ap¬ 
proximately. It is thus likely that the line widths 
are dependent on an interaction between the pri¬ 
mary and secondary. The variation could qualita¬ 
tively be explained as the result of a cavity in the 
primary wind produced by wind-wind collision. 
According to this model, the cavity created by the 
colliding winds would be of comparatively lower 
density, causing a reduction in blueward or red- 


ward emission at conjunctions, and in principle 
making the lines narrower at conjunctions than at 
quadratures of the binary phase. The spectra pre¬ 
sented in this paper are possibly consistent with 
this idea, although additional short-term variabil¬ 
ity of the line widths is suggested. 

One goal of the 2012 Chandra observing pro¬ 
gram of 6 Ori Aa was to allow observations of a 
massive star stellar wind as the short-period sec¬ 
ondary occulted different regions of emission for¬ 
mation on its journey around the primary star. We 
predict the reduction in the flux levels due to oc- 
cultation to be about 1-3% at most. Additional 
variability from other sources of greater magni¬ 
tude, as well as limited signal-to-noise in the data, 
make it impossible to identify occultation in our 
dataset at such a low percentage when there are 
clearly variations in the 10-15% range. In par¬ 
ticular, detailed analysis of Mg XII showed that 
flux, radial velocity, and FWHM vary both within 
a single orbit and within the dataset as a whole. 

The variability we see in the emission from 
6 Ori Aa is probably a composite of several ef¬ 
fects, including the long-term, greater than 9-day, 
photometric variability, binary orbit FWHM ef¬ 
fects, inter-orbit variability and intra-orbit vari¬ 
ability. It is likely that CIRs and/or pulsations 
play an important role in the variability. New long 
observations with the higher sensitivities offered 
by XMM-Newton would probably help resolve 
some of the photometric issues. Questions re¬ 
main concerning the source of the periods, phase- 
dependency of line profiles, various time scales 
of variability, and detailed modeling of the line 
width variability. Chandra observations at specific 
phases, such as conjunction and quadrature, and 
with a longer timeline, would be useful in verify¬ 
ing the model as well as parameterizing the vari¬ 
ability we have seen. Additional analysis of the 
UV DACs may clarify the sources of some of the 
components of the variability and in particular the 
rotation period of <5 Ori Aa. 
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Table 11 

Silicon Line Ratios and Derived Parameters 


phase 

MJD 

/// ratio 

G ratio 

77/77 e ratio 

Rfir/R* 

Ta MK 

ThHe 

MK 




Time Ordered 






.646 

56280.93 

0.52 ±0.26 

0.41 ±0.12 

0.25 ±0.09 

< 1.02 


8.23 

± 

0.79 

.734 

56281.38 

1.30 ±0.46 

0.83 ±0.17 

0.38 ±0.11 

1.46 ± 0.41 

> 4.50 

9.36 

± 

0.83 

.111 

56281.84 

0.73 ± 0.37 

1.47 ±0.72 

0.41 ±0.24 

1.10 ±0.09 

6.45 ±3.99 

9.27 

± 

1.66 

.082 

56283.49 

1.19 ±0.35 

0.94 ±0.19 

0.44 ±0.11 

1.32 ±0.28 

7.16 ±3.36 

9.73 

± 

0.74 

.170 

56283.97 

3.21 ± 1.14 

0.87 ±0.20 

0.20 ±0.09 

> 2.5 

10.20 ±6.18 

7.82 

± 

0.87 

.257 

56284.44 

1.99 ±0.96 

0.42 ± 0.09 

0.32 ± 0.08 

> 1.1 


8.93 

± 

0.70 

.475 

56285.79 


0.49 ± 0.09 

0.28 ±0.07 



8.60 

± 

0.57 

.562 

56286.28 

1.42 ± 0.40 

1.25 ±0.25 

0.32 ±0.11 

1.56 ± 0.42 

3.73 ±0.83 

8.90 

± 

0.94 

.649 

56286.77 

2.19 ±0.77 

0.94 ±0.19 

0.12 ±0.05 

> 1.5 

7.13 ±3.38 

6.99 

± 

0.59 

.955 

56288.40 

3.95 ±2.29 

0.57 ±0.14 

0.17 ±0.08 

> 1.7 

> 12 

7.49 

± 

0.83 

.042 

56288.89 

2.91 ± 1.14 

0.83 ±0.16 

0.46 ±0.10 

> 1.9 

11.11 ±6.39 

9.86 

± 

0.63 

.129 

56289.38 

0.99 ±0.22 

0.73 ±0.11 

0.25 ±0.07 

1.14 ±0.13 

> 8 

8.31 

± 

0.63 

Note. 

—Null entries imply unresolved ratio and/or parameter ranges. 






Table 12 

Magnesium Line Ratios and Derived Parameters 

phase 

MJD 

f /i ratio 

G ratio 

H/He ratio 

Rfir/R* 

Tg MK 

Thhs MK 





Time Ordered 




.646 

56280.93 

0.54 ±0.13 

0.94 ±0.16 

0.32 ±0.07 

2.30 ± 0.34 

4.35 ± 1.64 

5.71 ±0.36 

.734 

56281.38 

0.66 ±0.19 

0.93 ±0.18 

0.27 ±0.06 

2.61 ±0.49 

4.82 ± 2.12 

5.46 ±0.31 

.777 

56281.84 

0.51 ±0.15 

1.02 ±0.21 

0.48 ±0.10 

2.20 ±0.41 

3.86 ± 1.50 

6.33 ±0.33 

.082 

56283.49 

0.38 ±0.11 

0.96 ±0.18 

0.42 ± 0.08 

1.85 ±0.31 

4.20 ± 1.63 

6.14 ±0.32 

.170 

56283.97 

1.26 ±0.60 

0.62 ±0.16 

0.36 ±0.07 

4.55 ± 1.94 

36.40 ± 30.34 

5.89 ±0.36 

.257 

56284.44 

0.83 ± 0.23 

0.86 ±0.16 

0.44 ± 0.08 

3.10 ±0.62 

17.48 ± 14.41 

6.21 ±0.30 

.475 

56285.79 

1.05 ±0.53 

0.82 ±0.23 

0.38 ±0.08 

3.73 ± 1.49 

25.29 ± 22.37 

5.97 ±0.37 

.562 

56286.28 

0.42 ±0.13 

0.87 ±0.15 

0.37 ±0.07 

1.97 ±0.36 

16.95 ± 13.89 

5.97 ±0.35 

.649 

56286.77 

3.05 ± 1.67 

0.41 ±0.09 

0.26 ±0.05 

4.54 ± 0.00 

> 50 

5.43 ±0.24 

.955 

56288.40 

0.68 ±0.17 

1.19 ±0.21 

0.50 ±0.10 

2.71 ± 0.46 

2.68 ± 0.64 

6.42 ±0.28 

.042 

56288.89 

0.57 ±0.14 

1.11 ±0.18 

0.46 ± 0.09 

2.38 ±0.37 

3.02 ± 0.75 

6.27 ±0.30 

.129 

56289.38 

0.60 ±0.19 

0.82 ±0.17 

0.44 ± 0.08 

2.45 ± 0.50 

20.51 ± 17.27 

6.21 ±0.31 
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Table 13 

Neon Line Ratios and Derived Parameters 


phase 

MJD 

/// ratio 

G ratio 

H/He ratio 

Rf., 

ir/R* 

Tg MK 

ThHb 

MK 




Time Ordered 






.646 

56280.93 

0.23 ±0.10 

0.60 ±0.11 

0.82 ±0.12 

3.93 

±0.96 


3.66 ± 

0.17 

.734 

56281.38 

0.21 ±0.08 

0.80 ±0.13 

0.70 ±0.10 

3.66 

±0.73 

> 2.7 

3.49 ± 

0.14 

111 

56281.84 

0.19 ±0.06 

1.55 ±0.28 

1.01 ±0.17 

3.54 

±0.61 

1.05 ±0.17 

3.89 ± 

0.21 

.082 

56283.49 

0.17 ±0.05 

1.19 ±0.18 

1.07 ±0.14 

3.35 

±0.60 

1.63 ± 0.49 

3.97 ± 

0.16 

.170 

56283.97 

0.27 ±0.11 

0.56 ±0.11 

0.92 ±0.13 

4.24 

± 1.00 


3.80 ± 

0.18 

.257 

56284.44 

0.25 ±0.12 

0.40 ±0.10 

0.66 ±0.10 

4.02 

± 1.15 


3.44 ± 

0.14 

.475 

56285.79 

0.33 ±0.11 

0.52 ±0.09 

0.71 ±0.09 

4.86 

±0.91 


3.51 ± 

0.13 

.562 

56286.28 

0.14 ±0.08 

0.93 ±0.20 

0.87 ±0.14 

2.89 

±0.96 

> 1.6 

3.73 ± 

0.19 

.649 

56286.77 

0.40 ±0.10 

0.87 ±0.15 

0.78 ±0.12 

5.39 

±0.77 

> 2.1 

3.60 ± 

0.17 

.955 

56288.40 

0.21 ±0.06 

0.99 ±0.16 

0.75 ±0.11 

3.65 

±0.59 

2.57 ± 1.05 

3.57 ± 

0.15 

.042 

56288.89 

0.31 ± 0.08 

1.03 ±0.15 

0.85 ±0.11 

4.74 

±0.67 

2.26 ±0.82 

3.70 ± 

0.16 

.129 

56289.38 

0.38 ±0.13 

0.54 ±0.10 

0.83 ±0.11 

5.16 

±0.96 


3.67 ± 

0.15 

Note. 

—Null entries imply unresolved ratio and/or parameter ranges. 






Table 14 

Oxygen Line Ratios and Derived Parameters 

phase 

MJD 

f/i ratio 

G ratio 

H/He ratio 

Rfir/R* 

To MK 

ThHb 

MK 




Time Ordered 






.646 

56280.93 

0.12 ±0.08 

0.84 ±0.20 

1.05 ±0.19 

9.05 ± 3.40 

3.40 ± 1.90 

2.32 

± 

0.13 

.734 

56281.38 

0.06 ± 0.05 

2.50 ± 1.01 

2.55 ±0.94 

6.29 ±2.88 

0.00 ± 0.62 

2.98 

± 

0.36 

.777 

56281.84 

< 0.08 

0.96 ±0.27 

1.26 ±0.31 

< 6.9 

2.70 ± 1.63 

2.43 

± 

0.18 

.082 

56283.49 

< 0.02 

1.19 ±0.34 

1.73 ±0.39 

< 3.8 

1.59 ±0.99 

2.68 

± 

0.17 

.170 

56283.97 

< 0.02 

0.83 ± 0.25 

1.50 ±0.32 

< 4.0 

4.73 ± 3.32 

2.56 

± 

0.16 

.257 

56284.44 

< 0.02 

0.90 ±0.27 

1.43 ±0.31 

< 4.1 

3.31 ± 2.14 

2.52 

± 

0.16 

.475 

56285.79 

< 0.02 

0.74 ±0.23 

1.17 ±0.25 

< 3.9 

> 1.8 

2.39 

± 

0.15 

.562 

56286.28 

< 0.04 

0.68 ±0.18 

1.37 ±0.26 

< 5.3 

> 2.5 

2.50 

± 

0.14 

.649 

56286.77 

0.09 ± 0.08 

0.36 ±0.11 

0.79 ±0.14 

6.52 ±3.89 

> 77 

2.15 

± 

0.09 

.955 

56288.40 

< 0.08 

1.11 ±0.30 

1.66 ±0.35 

< 7.6 

1.78 ± 1.07 

2.64 

± 

0.16 

.042 

56288.89 

0.05 ± 0.05 

0.97 ±0.45 

1.60 ±0.53 

5.26 ±3.31 

> 0.35 

2.59 

± 

0.26 

.129 

56289.38 

< 0.04 

0.83 ± 0.22 

1.38 ±0.27 

< 7.2 

3.70 ± 2.22 

2.50 

± 

0.14 
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